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2.  Objectives 


The  objectives  are  to  create  and  evaluate  theoretical  and  experimental 
methods  for  active  control  of  nonlinear  aeroelastic  systems. 

3.  Status  of  effort 


Considerable  progress  has  been  made  toward  achieving  the  research 
objectives.  Both  mathematical  and  experimental  models  have  been 
completed,  and  excellent  agreement  has  been  achieved  between  results  from 
the  two  complementary  approaches.  The  uncontrolled  model  has  been  fully 
evaluated  and  initial  results  obtained  for  the  model  with  active  control  as 
well.  This  project  is  right  on  schedule. 

4.  Accomplishments/New  Findings 

The  successful  correlation  between  theory  and  experiment  has  resulted 
in  the  production  of  four  journal  papers  (see  below).  As  mentioned  above, 
the  models  are  now  being  used  to  develop  and  assess  various  control 
strategies  which  aim  to  increase  the  flutter  speed  of  this  aeroelastic  system. 
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the  department  of  Mechanical  Engineering  and  Materials  Science.  Part  of 
this  research  has  formed  the  basis  for  his  PhD  dissertation  which  was 
completed  this  past  summer. 
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Accepted  for  publication  in  the  Journal  of  Fluids  and 
Structures. 

3)  K.D.  Frampton,  R.L.  Clark  and  E.H.  Dowell  (1996)  "State- 
Space  Modeling  for  Aeroelastic  Panels  with  Linearized 
Potential  Flow  Aerodynamic  Loading".  Accepted  for 
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8.  New  discoveries,  inventions,  or  patent  disclosures 

Initial  indications  are  that  the  flutter  speed  of  this  type  of  airfoil  can 
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Active  Control  of  Panel  Flutter  with  Piezoelectric  Transducers 


Kenneth  D.  Frampton,*  Robert  L.  Clark, f  and  Earl  H.  Dowell^ 

Duke  University,  Durham,  North  Carolina  27708-0302 

This  article  investigates  the  active  control  of  panel  flutter  with  piezoelectric  transducers  and  including 
linearized  potential  flow  aerodynamics.  The  aerodynamic  modeling  is  accomplished  by  approximating  the 
aerodynamic  generalized  forces  with  infinite  impulse  response  filters.  These  filters  are  coupled  to  the  in 
vacuo  panel  dynamic  system  in  feedback,  thus,  creating  a  coupled,  aeroelastic  system.  The  panel  model 
is  developed  from  a  Rayleigh -Rltz  approach  and  includes  the  mass  and  stiffness  effects  of  a  piezoelectric 
transducer.  Acting  as  a  self-sensing  actuator,  the  piezoelectric  transducer  is  used  to  implement  direct  rate 
feedback  control.  Results  of  an  analytical  implementation  of  this  control  system  demonstrate  a  significant 
increase  in  the  flutter  boundaries. 


Nomenclature 

A,  B,  F,  G  =  state-space  matrices 

a  =  panel  length,  chord 

cu  =  speed  of  sound 

b  =  panel  width,  span 

Cp  =  piezoelectric  capacitance  matrix 

D  =  £7/(1  —  v^),  stiffness 

D,,,,  =  aerodynamic  influence  coefficient 

=  piezoelectric  constant 
77.™(0  =  aerodynamic  infiuence  function 

h  =  thickness 

ImniO  =  aerodynamic  influence  function 

i  =  output  current  vector 

K  =  stiffness  matrix 

M  =  Mach  number 

Mp  =  structural  mass  matrix 

M,  =  piezo  mass  matrix 

N  =  number  of  expansion  modes 

P(X)  =  approximate  z-transfer  function 

p(x,  y,  t)  =  aerodynamic  pressure 

QrO)  =  generalized  force 

9n(0  =  generalized  coordinate 

r„(i)  =  qn(,s)lh, 

=  aerodynamic  influence  coefficient 
Sp  =  piezo  area  moment  of  inertia 

s  =  tUJa 

t  =  time 

Up  =  freestream  velocity 

u  =  state-space  input  vector 

V  = 

Vmt,  =  piezo  saturation  voltage 

v„  =  input  voltage 

wCx,  y,  0  =  panel  displacement 

x,y,z  =  Cartesian  coordinates 

X  =  state  vector 

y  =  output  vector 

a,  -y  =  Fourier  transform  variables 

£  =  d^,v,pp„a^EpSplDph,{\  -  Vp) 
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0  =  electromechanical  coupling  matrix 

A  =  PpUW/D 

=  PcOlpsh. 

Fp  =  Pphp/p.h, 

p  =  density 

tr  =  Dp/D, 

T  =  integration  dummy  variable 

<l>(x,  y,  z,  t)  =  fluid  velocity  potential 

y)  =  modal  expansion  functions 

Subscripts 

a  =  aerodynamic 

m,  n  =  modal  indices 

p  =  piezoelectric 

s  =  structural 

Introduction 

HE  recent  emergence  of  technology  capable  of  perform¬ 
ing  the  necessary  computational  and  actuation  tasks  has 
led  to  the  application  of  active  control  to  many  problems  in 
dynamics.  However,  while  much  attention  has  been  focused  in 
the  literature  on  active  control  of  wing  flutter,  very  little  has 
been  done  in  the  area  of  panel  flutter  control.  The  phenomena 
of  panel  flutter  has  been  studied  for  many  years,  but  the  few 
available  references  that  apply  active  control  to  panel  flutter 
use  piston  theory  or  quasisteady  aerodynamics  as  the  working 
model. The  use  of  piston  theory  is  not  without  good  reason 
since  it  is  reasonably  accurate  at  higher  Mach  numbers  and 
modeling  the  effects  of  linearized  potential  flow  (which  is  nec¬ 
essary  for  Mach  numbers  below  1.5),  in  a  state-space  repre¬ 
sentation  suitable  for  control  investigations  is  difficult.  How¬ 
ever,  it  is  the  transonic  and  low  supersonic  range  that  is  usually 
most  critical  for  flutter.  For  this  reason,  the  inclusion  of  full 
potential  flow  aerodynamics  is  important  to  effective  control 
system  design. 

This  article  investigates  the  active  control  of  panel  flutter 
including  linearized  potential  flow  aerodynamics.  Develop¬ 
ment  of  the  aeroelectroelastic  model  consists  of  combining  the 
electroelastic  panel  model  developed  by  Hagood  et  al.^  with 
the  aeroelastic  panel  model  introduced  by  Frampton  et  al.‘  The 
electroelastic  panel  model  combines  the  dynamics  of  a  plate 
with  the  dynamics  of  piezoelectric  transducers  through  a  Ray¬ 
leigh -Ritz  formulation.  The  aeroelastic  panel  model  is  created 
by  coupling  a  set  of  aerodynamic  force  approximating  filters 
in  feedback  with  an  in  vacuo  panel  model.  The  set  of  filters 
is  constructed  from  the  aerodynamic  influence  functions  as  de¬ 
veloped  by  Dowell,^'®  using  Prony’s  method. 

The  control  scheme  used  in  this  investigation  is  direct  rate 
feedback  (DRFB).  This  is  implemented  by  using  a  piezoelec- 
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trie  transducer  simultaneously  as  a  sensor  and  actuator  as  de¬ 
scribed  by  Anderson  et  al.,’  Dosch  and  Inman,®  and  more  re¬ 
cently  Cole  and  Clark.’  The  piezoelectric  transducer  sensor 
signal,  which  is  proportional  to  the  strain  rate,  is  amplified  and 
fed  back  to  itself  as  a  control  signal. 

The  effectiveness  of  this  control  system  in  delaying  the  onset 
of  flutter  is  presented.  In  addition,  the  importance  of  consid¬ 
ering  full  potential  flow  aerodynamics  in  a  control  system  de¬ 
sign  is  demonstrated.  Also,  the  effect  of  piezoelectric  trans¬ 
ducer  saturation  is  addressed. 


Aeroelectroelasticity 
Panel  and  Piezoelectric  IVansducer  Dynamics 

The  panel  in  this  investigation  combines  the  aeroelastic 
modeling  introduced  by  Frampton  et  ah'*  with  the  electroelastic 
modeling  described  by  Hagood  et  al.®  This  method  employs  a 
Rayleigh-Ritz  formulation  to  discretize  the  coupled  equations 
of  motion  for  the  piezoelectric  sensor/actuator  and  the  aero¬ 
elastic  panel.  The  combination  of  the  two,  which  includes  the 
mass  and  stiffness  effects  of  the  panel  and  the  piezoelectric 
transducer,  is  referred  to  as  an  aeroelastic  piezostructure. 

This  method  assumes  a  separable  solution  using  the  in  vacuo 
panel  eigenfunctions  and  generalized  coordinates  of  the  form 

N 

w{x,  y,  t)  =  X  yMO  (1) 

n-I 

and  a  linear  voltage  distribution  across  the  piezoelectric  thick¬ 
ness.  The  resulting  actuation  equations  of  motion  with  no 
structural  damping  are 

[M,  +  +  [K,  +  K,]{q]  =  [0]{v}  -  pMliQ]  (2) 

and  the  rate  sensor  equations  are 

+  [Cp]{F}  =  {^}  (3)  and  a  finiteness  or  radiation  condition  as  z  approaches  infinity. 

The  solution  can  be  obtained  by  taking  a  Laplace  transform 
where  with  respect  to  time  and  a  double  Fourier  transform  with  re¬ 

spect  to  the  X  and  y  spatial  coordinates.  The  transformation  is 


II 

(4) 

applied  to  Eqs.  (10),  (11),  and  Bernoulli’s  equation; 

{Q)  =  {Q>  2:  •••  Qnf 

{v)  =  {v,  V2 

(5) 

(6) 

,  ,, 

(12) 

The  aerodynamic  generalized  forces  Q„  are  a  function  of  the 
aerodynamic  pressure 

Qn{t)  =  f  f  — ^„{x,  y)  dx  dy  (7) 

Jo  Jo  PaUc 

while  incorporating  Eq.  (1).  Details  of  the  solution  process  can 
be  found  in  Refs.  5  and  6. 

The  solution  yields  the  generalized  aerodynamic  forces  on 
the  panel  that  are  given  here  as 

N 

Definitions  of  the  previous  matrices  can  be  found  in 
The  coordinate  system  for  the  panel  is  shown  in  Fig. 

Ref.  3. 

1. 

Qnit)  =  2  2™ 

mal 

(13a) 

Full  Potential  Flow  Generalized  Forces 

where 

Consider  the  case  of  a  simply  supported  panel  with  the  fol¬ 
lowing  eigenfunctions  (only  one  mode  in  the  y  direction  is 
used  for  simplicity): 

2™(0  =  +  qM)D^  +  1  qJ,r)H^{t  -  t)  dr 

y)  =  sin[(n7r/£i)x]sin[(7r/I»)y] 

(8) 

Jo 

r 

In  this  case  the  panel  equation  of  motion  [Eq.  (2)] 
expressed  nondimensionally  as 

can  be 

+  1  -  t)  dr 

(13b) 

[M,  +  +  {r} 

(13c) 

=  8  ^  [0]{v}  -  4/Xa{e} 

A 

(9) 

~  A/t/.  [  1 

(13d) 

Fig.  1  Panel  coordinate  system. 


where  the  hat  [''],  indicates  that  all  dimensional  terms  have 
been  pulled  out  of  the  matrix  and  included  in  the  nondimen- 
sional  coefficients. 

The  generalized  forces  are  found  by  solving  the  partial 
differential  equation  describing  the  velocity  potential  in  an  in- 
viscid,  irrotational  fluid  moving  parallel  to  the  x  axis 

_3  1  r  d^d>  d^<f>  , 

^  =0  (10) 
aL  Ldt  dxdt  dx  W 


subject  to  the  boundary  conditions  for  a  panel  embedded  in  an 
infinite  baffle 


dw  ,  ,,  dw 

- 1-  t/a  —  on  the  panel 

dt  dx 

0  off  the  panel 
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X  +  "y^]  da  dy  (13e) 

■  -S3?  £  £, 

X  J,[a^t\^a^  +  y^]  da  d-y  (13f) 

G™  =  I  I  'I^„(;c,  dy  dx 

X  [  [  'I'„(x,  dy  dx  (13g) 

Jo  Jo 

The  integrals  in  Eqs.  (13c),  (13d),  and  (13g)  can  be  per¬ 
formed  analytically  for  most  panel  eigenfunctions.  Equations 
(13e)  and  (13f),  which  define  the  aerodynamic  influence  func¬ 
tions,  must  be  evaluated  numerically. 


Fig.  2  State-space  schematic  of  the  complete  aeroelastic  system. 

approximations  to  the  influence  functions  can  be  readily  trans¬ 
formed  into  an  /th  order  infinite  impulse  response  (HR)  filter 
of  the  form 


Piston  Theory  Generalized  Forces 
Piston  theory  is  presented  here  for  comparison  with  the  full 
potential  flow  solution  developed  previously.  Piston  theory  as¬ 
sumes  that  the  pressure  acting  on  Ae  panel  is  equivalent  to  the 
pressure  acting  on  a  piston  in  a  tube; 


P  =  P«a- 


dw  .  dw 


(14) 


The  total  piston  velocity  includes  a  convection  term  U„(,dwl 
dx),  as  well  as  a  direct  velocity  term  dw/dt.  Here  the  plate  is 
the  equivalent  piston  and  the  tube  is  perpendicular  to  the  plate. 

Combining  Eqs.  (1)  and  (14)  and  inserting  into  Eq.  (7) 
yields  the  piston  theory  generalized  aerodynamic  forces 


where 


N 


QM 

=22™ 

■■••1 

(15a) 

Q^d)  =  qMS,„„  +  q„(f)D™, 

(15b) 

(15c) 

■o 

dy  dx 

)  Jo 

(15d) 

Note  that  this  result  is  identical  to  that  obtained  from  full 
potential  flow  analysis  if  the  influence  functions  and 

/™(r)]  are  ignored. 

Approximate  Generalized  Forces 

Since  the  generalized  aerodynamic  forces  represented  in  Eq. 
( 1 3b)  are  a  function  of  the  panel  generalized  coordinates,  they 
can  be  viewed  as  dynamic  feedback  on  the  panel.  These  feed¬ 
back  dynamics  are  characterized  by  the  influence  functions 
//„„,(0  and  /„„,(/),  which  represent  aerodynamic  impulse  re¬ 
sponses,  and  the  influence  coefficients  S™  and  D™,,  which  are 
instantaneous  or  feed-through  dynamics. 

Since  there  is  no  closed-form  solution  for  the  generalized 
forces  some  approximation  must  be  made.  The  approach  sug¬ 
gested  here  is  to  construct  digital  filters  that  approximate  these 
dynamics.  A  reduced  order  recursive  filter  may  be  constructed 
by  applying  Prony’s  method  to  the  aerodynamic  influence 
functions.'""  Prony’s  method  approximates  an  impulse  re¬ 
sponse  (in  this  case  the  influence  functions)  with  a  set  of  ex¬ 
ponential  functions  in  a  least-squares  sense.  These  exponential 


(16) 


As  an  example,  suppose  the  values  of  the  influence  function 
are  used  to  solve  for  the  filter  coefficients  in  Eq.  (16). 
When  provided  with  an  input  signal  that  represents  the  panel 
generalized  coordinate  q„{iT),  this  filter  would  output  an  ap¬ 
proximate  convolution  of  qJiJT)  with  Similarly,  a  fil¬ 

ter  could  be  constructed  from  I„„(iT)  and  provided  with  a  sig¬ 
nal  representing  Summing  the  outputs  of  the  two  filters 

would  yield  an  approximation  to  the  dynamic  portion  of  Eq. 
(13b).  Details  of  this  process  may  be  found  in  Ref.  4. 

An  approximate  filter  can  be  found,  as  shown  previously, 
for  each  aerodynamic  influence  function  and  /m,(f).  A 

complete  set  of  these  filters  can  then  be  summed  according  to 
Eqs.  (13a)  and  (13b).  This  would  yield  a  dynamic  system  that 
has  the  generalized  coordinates  q„{t)  and  ^„(t)  as  the  inputs 
and  the  generalized  aerodynamic  forces  Q„(t)  as  the  outputs. 
This  system  is  then  coupled  as  feedback  on  the  panel  resulting 
in  a  complete  aeroelastic  system.  This  is  shown  schematically 
(in  state-space  notation)  in  Fig.  2. 

State-Space  Formulation 

The  actuation  equations  of  motion  for  the  panel  described 
in  Eq.  (2)  can  be  expressed  in  state  variable  form  as 


where 


X,  =  A,x,  +  B,u, 
y.  =  F,x, 


(17a) 

(17b) 


S.= 


PaUl[M,  +  A/,]-’  [Af, 


F,  = 


1 

& 


3-  =  ^ 

l«oJ 

1  (17c) 

Af,]  ( 

;] 

(17d) 

0 

+  A/,] 

-10 

(17e) 

(17f) 

The  approximate  filters  of  Eq.  (16)  can  be  converted  to  state- 
space  form.'^  Each  approximate  filter  associated  with  the  in¬ 
fluence  function  H„Jj)  and  /„„,(0  has  a  state-space  represen¬ 
tation  with  matrices  {A^^,  F„^,  GhJ  and 

respectively.  The  complete  set  of  state-space  repre- 
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sentations  can  be  assembled  to  form  the  approximate  aerody¬ 
namic  system  as  follows: 


where 


=  AaXa  + 
ya  =  FcXa  -t-  G^a 


ya=[QM  Q,{t) 


A„  = 


s«  = 


F.  = 


G.= 


92 

92 

•••  9. 

9»1 

T 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A,^_ 

0 

0 

0 

1 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

J 

F,n 

F.S 

Fh„ 

F,^ 

Fh^ 

.Fh., 

F,.. 

Fh^ 

F, 

•mttrni 

'cl 

g;„ 

CL 

CL 

CL 

C';,^ 

CL 

CL 

CL 

CL 

CL. 

(18a) 

(18b) 

(18c) 

(18d) 


(18e) 


(18f) 


(18g) 


(18h) 


and  the  influence  coefficients  are  incorporated  as  follows: 


Gfj^  —  Gff^  +  Sn„  (18i) 

CL  =  +  D™  (18j) 

The  aerodynamic  states  x„  are  of  mathematical  construct  and 
have  no  physical  significance.  The  panel  and  aerodynamic 
state-space  models  can  then  be  coupled  together  as  shown  in 
Fig.  2. 

The  accuracy  of  this  model  can  be  demonstrated  by  com¬ 
paring  the  flutter  boundaries  with  those  reported  in  the  litera- 
ture.’  ®  '^  '"  The  flutter  boundaries  are  determined  by  increasing 
the  nondimensional  dynamic  pressure  A  while  holding  all  other 
parameters  fixed.  Flutter  occurs  when  one  of  the  system  ei¬ 
genvalues  moves  into  the  right-half-plane. 

A  comparison  of  flutter  dynamic  pressure  vs  Mach  number 
is  shown  in  Fig.  3.  Results  are  illustrated  for  a  panel  with  mass 
ratios  of  =  0. 1  and  /ij,  =  0.5,  cr  =  1 .2,  e  =  130,  and  dimen¬ 
sion  ratio  alb  =  1,  using  a  four-mode  panel  model.  The  figure 
shows  the  flutter  boundaries  for  the  plate  with  and  without  a 
piezoelectric  transducer  for  both  full  potential  flow  and  piston 
theory  aerodynamics.  Note  that  the  presence  of  the  piezoelec¬ 
tric  transducer  slightly  increases  the  flutter  boundary.  This  is 
because  of  the  increased  panel  mass  with  the  transducer.  Also 
note  that  the  two  aerodynamic  theories  approach  each  other 
asymptotically  as  the  Mach  number  increases.  However,  for 
Mach  numbers  below  1.5  the  theory  incorporating  linearized 
potential  flow  aerodynamics  predicts  flutter  boundaries  well 


1.0  1.2  1.4  1.6  1.8  2.0 

Mach  Number 

Fig.  3  kfysM  for  a  panel  aspect  ratio  of  alb  =  1. 

below  that  predicted  by  piston  theory.  These  results  are  con¬ 
sistent  with  those  obtained  by  Dowell. A  more  detailed  dis¬ 
cussion  of  the  model  validity  can  be  found  in  Ref.  4. 

Collocated  Direct  Rate  Feedback  Control 

Collocated  direct  rate  feedback  control  is  a  well-documented 
form  of  control  which,  in  the  absence  of  sensor/actuator  dy¬ 
namics  and  rigid  body  modes,  has  guaranteed  stability  for  any 
positive  feedback  gain.’’ 

Collocation  is  achieved  by  using  a  piezoceramic  transduc¬ 
tion  device  as  a  sensoriactuator,^  i.e.,  simultaneously  as  a  sen¬ 
sor  and  actuator.’'^'"*  By  applying  the  appropriate  circuitry  to 
the  piezoelectric  transducer,  the  portion  of  the  piezoelectric 
output  that  is  proportional  to  the  structural  velocity  can  be 
isolated.  This  is  equivalent  to  rendering  Eq.  (3)  as 

[©]^{4)  =  {1}  (19) 

Direct  rate  feedback  control  consists  of  feeding  back  a  signal 
proportional  to  the  sensor  current  i  from  Eq.  (19)  as  the  voltage 
input  V  to  Eq.  (2). 


Results 

The  effectiveness  of  collocated  direct  rate  feedback  control 
is  demonstrated  by  applying  the  control  to  the  aeroelastic  panel 
model  and  finding  the  flutter  instability  boundaries.  All  results 
in  this  section  are  for  the  following  parameters:  /x,  =  0.5, 

=  0.1,  cr  =  1.2,  c  =  130,  and  alb  =  1.0.  These  parameters 
correspond  to  typical  values  for  G-1 195-type  piezoceramic  and 
a  steel  panel.  Previous  investigations  have  shown  that  mass 
ratio  has  relatively  little  effect  on  the  flutter  boundaries,''* 
therefore,  the  mass  ratios  /x,  and  /x„  are  held  constant. 

Selection  of  the  piezoelectric  transducer  location  is  very  im¬ 
portant  for  a  successful  control  system  design.  Formal  opti¬ 
mization  of  transducer  location  was  not  performed,  however, 
investigation  into  the  aeroelastic  mode  shapes  suggests  a  good 
location.  It  is  well  known  that  piezoelectric  transducers  are 
most  effective  when  located  in  regions  of  high  strain  (i.e.,  high 
panel  curvature).  Since  the  goal  is  to  control  flutter,  the  flutter 
mode  shapes  must  be  considered  when  looking  for  advanta¬ 
geous  locations.  This  approach  was  also  considered  by  Scott 
and  Weisshaar.' 

The  normalized  flutter  mode  shapes  in  the  chordwise  direc¬ 
tion  are  shown  in  Fig.  4  for  various  Mach  numbers.  Note  the 
strong  influence  of  the  1st  and  2nd  in  vacuo  modes  on  the 
flutter  mode.  The  region  of  maximum  strain  at  each  Mach 
number  occurs  very  nearx/a  =  0.7.  For  this  reason  a  transducer 
coverage  of  0.6  <  xia  <  0.8  and  0.4  <  ylb  <  0.6  was  selected. 
This  chordwise  location  also  corresponds  to  the  maximum  de¬ 
lay  in  the  migration  of  open  loop  zeros  into  the  right-half- 
plane.  The  importance  of  the  zero  migration  delay  will  be  dis¬ 
cussed  in  more  detail  later. 
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Results  from  this  analysis  are  shown  in  Fig.  5  for  the  full 
potential  flow  model  and  several  feedback  gains.  Figure  5 
demonstrates  that,  except  near  M  =  1.0,  the  use  of  direct  rate 
feedback  control  is  equally  effective  at  all  Mach  numbers  in 
terms  of  a  percentage  increase  in  flutter  dynamic  pressure. 
Note  in  Fig.  5  that  there  is  a  maximum  flutter  boundary  that 
direct  rate  feedback  can  achieve.  This  is  demonstrated  by  the 
coincidence  of  the  flutter  boundary  for  feedback  gains  of  0. 1 
and  0.2. 

This  behavior  is  better  illustrated  in  Fig.  6,  which  shows  the 
change  in  flutter  nondimensional  dynamic  pressure  with  in¬ 
creasing  feedback  gain  for  three  Mach  numbers.  Also  included 
in  Fig.  6  are  similar  results  for  a  piston  theory  aerodynamic 
model.  It  is  apparent  in  Fig.  6  that  the  piston  theory  model 
overpredicts  the  control  effectiveness,  particularly  at  lower 
Mach  numbers,  compared  to  the  full  potential  flow  model. 
Also  note  that  the  difference  between  the  full  potential  flow 
model  and  the  piston  theory  model  increases  with  increasing 
gain. 

An  important  result  demonstrated  in  Fig.  6  is  that,  at  a  fixed 
Mach  number,  there  is  a  maximum  flutter  dynamic  pressure 
that  direct  rate  feedback  control  can  achieve.  This  is  evident 
in  the  asymptotic  behavior  of  the  flutter  boundaries  with  in¬ 
creasing  gain  shown  in  Fig.  6.  In  other  words,  for  a  fixed  Mach 
number,  there  is  a  limiting  nondimensional  dynamic  pressure 
A,  beyond  which  direct  rate  feedback  control  cannot  stabilize 
the  system. 

This  Umitation  is  because  of  the  mechanism  by  which  direct 
rate  feedback  control  stabilizes  the  system.  Stabilization  is  ac¬ 
complished  by  moving  the  open-loop  poles  toward  the  open- 
loop  zeros.  'nUs  is  demonstrated  in  Fig.  7,  which  shows  the 
closed-loop  root  locus  for  a  four-mode  panel  (aerodynamic 


Plate  Chord,  xl  a 

Fig.  4  Flutter  mode  shapes  for  various  Mach  numbers. 


poles  are  omitted  for  clarity).  Results  in  Fig.  7  are  for  a  Mach 
number  of  2.0  and  dynamic  pressure  A  =  1300,  which  is  greater 
than  the  open-loop  flutter  dynamic  pressure,  but  less  than  the 
limiting  dynamic  pressure.  Note  in  Fig.  7  that  as  the  feedback 
gain  is  increased  the  unstable  poles  are  drawn  toward  the  open- 
loop  zeros  that  are  located  in  the  left-half-plane.  Therefore, 
given  enough  feedback  gain,  the  system  can  be  stabilized. 

Unfortunately,  it  is  not  always  the  case  that  the  open-loop 
zeros,  toward  which  the  unstable  poles  are  drawn,  lie  in  the 
left-half-plane.  The  location  of  these  open-loop  zeros  is  a  func¬ 
tion  of  M,  mass  ratio  the  sensoriactuator  location,  and  most 
strongly  a  function  of  the  dynamic  pressure  A.  This  depen¬ 
dence  is  demonstrated  in  Fig.  8,  which  shows  the  locus  of  the 
system  zeros  as  a  function  Of  A.  This  plot  is  for  a  Mach  number 
of  2.0.  Note  the  migration  path  of  the  zeros.  At  a  dynamic 
pressure  of  A  =  1593,  a  pair  of  zeros  moves  into  the  right-half¬ 
plane.  Piston  theory  pr^icts  that  the  zero  will  move  into  the 
right-half-plane  at  a  dynamic  pressure  of  A  =  1903.  This  ac¬ 
counts  for  the  large  difference  in  maximum  flutter  dynamic 
pressure  depicted  in  Fig.  6. 

A  root  locus  similar  to  that  in  Fig.  7,  but  with  A  =  1630,  is 
shown  in  Fig.  9.  Note  how  the  pair  of  unstable  poles  is  moved 
toward  a  pair  of  zeros  in  the  right-half-plane.  Therefore,  no 
amount  of  feedback  gain  could  stabilize  the  system.  The  point 
at  which  the  open-loop  zeros  moves  into  the  right-half-plane 
is  the  closed-loop  flutter  boundary.  Beyond  this  point,  even 
with  direct  rate  feedback  control,  the  panel  will  flutter. 


Feedback  Gain 


Fig.  6  Flutter  boundaries  with  varying  feedback  gain. 


Fig.  7  Root  locus  of  the  full  potential  flow  model  (A  =  1300)  with 
feedback  gain  varying  from  0  (corresponding  to  the  open-loop 
pole,  x)  to  »  (corresponding  to  the  open-loop  zero,  o). 


Fig.  8  Zero  locus  of  the  full  potential  flow  model  with  nondi- 
mensional  dynamic  pressure  varying  from  A  =  0  to  A  =  1600 
(o  =  zero). 


Fig.  9  Root  locus  of  the  full  potential  flow  model  (A  =  1630)  with 
feedback  gain  varying  from  0  (corresponding  to  the  open-loop 
pole,  X)  to  (corresponding  to  the  open-loop  zero,  o). 

Although  this  demonstrates  a  limit  to  the  amount  that  direct 
rate  feedback  control  can  increase  the  flutter  boundary,  other 
feedback  control  techniques  could  potentially  increase  the 
boundary  further.  Since  the  migration  of  the  open-loop  zeros 
is,  in  part,  a  function  of  the  location  of  the  sensoriactuator,  its 
size  and  location  on  the  aeroelastic  panel  could  be  formally 
optimized  to  further  delay  the  zero  migration  into  the  right- 
half-plane. 

A  second  limitation  in  the  effectiveness  of  DRFB  flutter  con¬ 
trol  is  the  saturation  of  the  piezoelectric  transducer.  The  ap¬ 
plied  control  signal  voltage  is  proportional  to  the  strain  rate 
and  is  therefore  a  function  of  the  panel  disturbance.  In  a  re¬ 
alistic  situation  the  disturbance  would  result  from  externally 
applied  pressures  such  as  turbulent  boundary-layer  (TBL)  pres¬ 
sures.  For  simplicity  in  this  investigation  an  initial  displace¬ 
ment  in  the  first  mode  of  r,  =  qjh,  =  1.0  was  used  as  the 
disturbance.  The  time  history  of  the  control  signal  was  cal¬ 
culated  for  various  Mach  numbers  and  the  dynamic  pressure 
at  which  the  maximum  allowable  piezoelectric  voltage  was 
reached  was  noted. 

Results  for  transducer  saturation  are  shown  in  Fig.  6.  Note 
that  for  M  =  1.2  the  saturation  curve  corresponds  with  the 
potential  flow  flutter  curve.  Not  until  higher  Mach  numbers 
does  the  piezoelectric  saturation  significantly  limit  the  ability 
to  delay  flutter  onset.  Considering  piezoelectric  saturation  ef¬ 


Frequency,  k  =  m[p^h^a"  /  D^Y' 

Fig.  10  Open-  and  closed-loop  transfer  functions. 

fects  the  flutter  dynamic  pressure  is  increased  by  55%  at  Af  = 
1.2,  28%  at  A/  =  1.5,  and  26%  at  M  =  2.0,  as  shown  in 
Fig.  6. 

Also  note  that  while  the  application  of  collocated  direct  rate 
feedback  control  is  effective  for  panel  flutter  control  it  is  also 
effective  for  disturbance  rejection.  This  is  demonstrated  in  Fig. 
10,  which  shows  the  open-  and  closed-loop  transfer  functions 
for  three  Mach  numbers.  In  each  case  the  dynamic  pressure  is 
just  below  the  corresponding  flutter  dynamic  pressure  for  that 
Mach  number.  The  effect  of  direct  rate  feedback  control  in  this 
case  is  to  add  damping  to  the  near  flutter  poles,  thus  reducing 
the  spectral  peak. 

Conclusions 

The  active  control  of  panel  flutter  has  been  investigated  in¬ 
cluding  full  linearized  potential  flow  aerodynamic  theory  ap¬ 
propriate  for  the  full  transonic  and  supersonic  Mach  number 
range.  Direct  rate  feedback  control  from  a  sensoriactuator  pi¬ 
ezoelectric  transducer  was  modeled  and  its  effectiveness  in 
controlling  flutter  instability  was  demonstrated.  It  was  shown 
that  this  type  of  control  is  capable  of  significantly  increasing 
the  flutter  nondimensional  dynamic  pressure.  It  was  also 
shown  that  this  type  of  control  is  limited  by  the  migration  of 
the  open-loop  zeros  into  the  right  half-plane,  thus,  beyond  a 
certain  dynamic  pressure,  preventing  the  system  from  being 
stabilized.  This  zero  migration  to  the  right  half-plane  defines 
a  new  flutter  boundary  for  the  closed-loop  aeroelastic  panel. 
A  second  limitation  was  demonstrated  in  the  saturation  of  the 
piezoelectric  transducer. 
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This  article  presents  a  new  method  for  the  state-space  modeling  of  aeroelastic  panels  subject  to  line¬ 
arized  potential  flow  aerodynamic  loading.  This  is  accomplished  by  approximating  the  aerodynamic  gen¬ 
eralized  forces  on  the  panel  with  discrete,  infinite  impulse  response  (UR)  filters.  The  UR  filters  are  created 
using  Prony’s  method  for  a  least-squares-fit  to  the  aerodynamic  influence  functions.  These  fiiters  are 
coupied  to  the  in  vacuo  panel  dynamic  system  as  feedback,  creating  a  coupled,  aeroelastic  system.  The 
accuracy  of  the  model  is  established  by  comparing  the  panel  flutter  boundaries  of  the  approximate  system 
with  those  found  in  past  studies. 
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Nomenclature 

=  state-space  matrices 
=  panel  length,  chord 
=  filter  coefficients 
=  speed  of  sound 
=  panel  width,  span 
=  filter  coefficients 
=  £//(!  —  i/^),  panel  stiffness 
=  aerodynamic  influence  coefficient 
=  aerodynamic  influence  function 
=  panel  thickness 

=  approximate  infinite  impulse  response  filter 
order 

=  aerodynamic  influence  function 
=  Bessel’s  function 
=  Mach  number 
=  modal  mass 

=  number  of  expansion  modes 
=  exact  z  transfer  function 
=  approximate  z  transfer 
=  approximate  Laplace  transfer  function 
=  aerodynamic  pressure 
=  generalized  force 
=  generalized  coordinate 
=  finite  impulse  response  filter  order 
=  qn{t)/h 

=  aerodynamic  influence  coefficient 
=  tUJa 

=  sample  period 
=  time 

=  freestream  velocity 
=  state-space  input  vector 
=  panel  displacement 
=  Cartesian  coordinates 
=  state  vector 
=  output  vector 
=  Fourier  transform  variables 
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K 

=  Laplace  domain  variable 

=  puy/D 

A 

M 

=  PcO/Pp 

Pc 

=  density  of  fluid 

Pp 

=  panel  mass  per  area 

T 

=  integration  dummy  variable 

(f)(x,  >',  z,  t) 

=  fluid  velocity  potential 

y) 

=  modal  expansion  functions 

=  modal  natural  frequency 

Subscripts 

a 

=  aerodynamic 

m,  n 

=  modal  indices 

P 

=  panel 

Introduction 

PANEL  flutter  is  the  self-excited  dynamic  instability  of 
plate-like  structures  exposed  to  fluid  flow.  This  aeroelastic 
phenomenon  has  received  much  attention  in  the  past  and  is 
very  important  toward  the  successful  development  of  super¬ 
sonic  aircraft  such  as  the  National  Aerospace  Plane,  the  High 
Speed  Civil  Transport,  and  the  Advanced  Tactical  fighter,  to 
name  a  few.  Dugundji'  published  an  excellent  paper  on  early 
investigations  into  linear  panel  flutter,  whereas  Gray  and  MeF 
gave  a  review  of  more  recent  efforts  including  nonlinear  ef¬ 
fects.  However,  the  vast  majority  of  work  has  been  performed 
assuming  piston  theory  as  the  aerodynamic  model.  The  use  of 
piston  theory  is  not  without  good  reason,  since  it  is  reasonably 
accurate  at  higher  Mach  numbers  and  modeling  the  effects  of 
linearized  potential  flow,  which  is  necessary  for  lower  Mach 
numbers  (M  <  1.5),  is  challenging.  Since  low  supersonic  flight 
speeds  are  utilized  by  many  aircraft,  it  is  important  to  consider 
the  effects  of  full  potential  flow  aerodynamic  loading  in  panel 
flutter  investigations. 

The  few  references  that  include  linearized  potential  flow  aer¬ 
odynamics  focus  on  the  determination  of  flutter  boundaries, 
and  subsequently,  the  limit  cycle  oscillations  in  the  postflutter 
(nonlinear)  regime.  These  models,  however,  are  not  conducive 
to  the  application  of  modem  control  theory.  A  method  of  build¬ 
ing  a  state-space  model  of  the  coupled  aeroelastic  system 
would  permit  the  determination  of  flutter  boundaries  as  well 
as  provide  a  means  of  applying  modem  control  theory  to  the 
problem.  While  this  has  already  been  accomplished  for  typical 
wing  sections,^"'  it  has  yet  to  be  done  for  panels. 

This  article  presents  a  method  for  constructing  a  state-space 
model  of  the  coupled  aeroelastic  panel  including  linearized 
potential  flow  aerodynamics.  This  is  accomplished  by  approx¬ 
imating  the  aerodynamic  generalized  forces  on  a  panel,  as  de¬ 
veloped  by  Dowell,^'*  with  discrete,  infinite  impulse  response 
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(UR) '  filters.  These  filters  are  coupled  to  the  in  vacuo  panel 
dynamic  system  as  feedback,  creating  a  coupled,  aeroelastic 
system.  The  accuracy  of  the  model  is  established  by  comparing 
the  panel  flutter  boundaries  of  the  approximate  system  with 
those  found  in  the  literature.*'’ 

Panel  Dynamics  and  Generalized  Forces 

The  partial  differential  equation  (PDE)  of  motion  describing 
a  thin,  uniform  panel  is‘° 

r  aV(x,  y,  t)  y,  t)  d*w(x,  y,  t) 

[  dx^  dx^dy^  dy* 

d^w(x,  y,  t) 

+  Pp  g/  +  Pix,  y,  0  =  0  (1) 

The  coordinate  system  for  the  panel  is  shown  in  Fig.  1. 

A  separable  expansion  solution  is  assumed  using  the  in 
vacuo  orthogonal  panel  eigenfunctions  and  generalized  coordi¬ 
nates; 


w(jt,  y,  t)  =  2  >’)?"(') 


Substituting  Eq.  (2)  into  Eq.  (1),  multiplying  by  an  arbitrary 
expansion  function  y),  and  integrating  over  the  domain 

yields  a  set  of  ordinary  differential  equations  of  the  form: 

+  p^UlQM  =  0  (3) 


=  |  I  dx  dy 

_  f 

-Jo  },[dx^  dx^ 

( a^'p„  a’^pA 


ax^  ay^  ax^  ay^ 


-(-  2(1  -  i/") 


a^T'„  a’^„ 

dxdy  dxdy 


a^a^' 

ay^  ay'  _ 


and  the  generalized  forces  are  a  function  of  the  aerodynamic 
pressure 


r  fpU.  y. 

.0  Jo 


■  T'„(x,  y)  dx  dy 


Full  Potential  Flow  Generalized  Forces 

Consider  the  case  of  a  simply  supported  panel  with  the  fol¬ 
lowing  eigenfunctions  (only  one  mode  in  the  y  direction  for 
simplicity  such  that  n  =  1,  2,  . . . ,  N): 

''I'n(^.  y)  =  sin[(n7r/a)x]sin[(7r/f>)y]  (7) 

In  this  case,  the  panel  equation  of  motion,  Eq.  (3),  can  be 
expressed  nondimensionally  as 

r„(s)  +  +  (a/bffr„(s)  +  4^lQ„(s)  =  0  (8) 

The  generalized  forces  Q„U)  are  found  by  solving  the  PDE 
describing  the  velocity  potential  in  an  inviscid,  irrotational 
fluid  moving  parallel  to  the  x  axis 


Fig.  1  Panel  coordinate  system. 

subject  to  the  boundary  conditions  for  a  panel  embedded  in  an 
infinite  baffle 


dw 

—  +  Ua  —  on  the  panel 
dt  dx 


U-o  0 


off  the  panel 


and  a  finiteness  or  radiation  condition  as  z  approaches  infinity. 

The  solution  can  be  obtained  by  taking  a  Laplace  transform 
with  respect  to  time  and  a  double  Fourier  transform  with  re¬ 
spect  to  the  X  and  y  spatial  coordinates.  The  transformation  is 
applied  to  Eqs.  (9)  and  (10)  and  Bernoulli’s  equation, 


(a^ 


while  incorporating  Eq.  (2).  Details  of  the  solution  are  dis¬ 
cussed  by  Dowell.** 

The  solution  yields  the  generalized  aerodynamic  forces  on 
the  panel,  which  are  given  here  as 

N 

Qn(t)  =  2  (12a) 

where  Q,^{t)  is  the  force  on  the  nth  panel  mode  because  of 
the  motion  of  the  mth  panel  mode,  such  that 

j2™(0  =  qJ.t)S^  +  4J.i)D^  +  j  q„(T)H^(t  -  t)  dr 
+  j  ‘}n,('r)/^(t  -  t)  dr  (12b) 

1  r°  r*  d'^ 

X  da  dy  (12e) 

X  Jt(a.ty/a^  +  y‘)  da  dy  (12f  ) 

X  f  f  ^„(x,  y)e"“^^>  dy  dx  (12g) 

Jo  Jo 


X  y,(a„t\/a'  +  •/)  da  dy 


/™(0  = 


G™Va'  -F 


X  -F  ■/)  da  dy 


y)e'''“"’»'  dy  dx 


yy^-^*yy>  jy  dx 
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The  integrals  in  Eqs.  (12c),  (12d),  and  (12g)  can  be  per¬ 
formed  analytically  for  most  panel  eigenfunctions.  Equations 
(12e)  and  (12f),  which  define  the  aerodynamic  influence  func¬ 
tions,  must  be  evaluated  numerically. 

Piston  Theory  Generalized  Forces 

Piston  theory  is  presented  here  for  comparison  with  the  full 
p>otential  flow  solution  developed  previously.  Piston  theory  as¬ 
sumes  that  the  pressure  acting  on  the  panel  is  equivalent  to  the 
pressure  acting  on  a  piston  in  a  tube: 

Tlie  total  piston  velocity  includes  a  convection  term  U„(dwl 
dx),  as  well  as  a  direct  velocity  term  dw/dt.  Here  the  plate  is 
the  equivalent  piston  and  the  tube  is  perpendicular  to  the  plate. 

Combining  Eqs.  (2)  and  (13)  and  inserting  into  Eq.  (6) 
yields  the  piston  theory  generalized  aerodynamic  forces 

N 

QM)  =  2  QUO  (14a) 

fR«l 

where 


represented  by  Eq.  (12b),  which  has  two  impulse  responses  of 
HUiT)  and  /™(jT).  Note  that  HUiT)  and  /™(jT)  are  discrete 
time  approximations  of  Eqs.  (12e)  and  (12f ). 

Equation  (15)  can  be  used  to  represent  the  aerodynamics  of 
Eq.  (12b)  by  substituting  a  numerically  calculated  influence 
function  [HUiT)  or  lUiT)]  for  h{iT)  in  Eq.  (15).  However, 
constructing  a  state-space  model  of  an  FIR  filter  requires  as 
many  states  as  filter  coefficients.  Therefore,  any  reasonable 
discrete  temporal  resolution  of  the  influence  functions  would 
result  in  a  very  high-order  state-space  model. 

Fortunately,  a  reduced  order  recursive  filter  may  be  con¬ 
structed  by  applying  Prony’s  method  to  the  influence  func¬ 
tions.'*'^  Prony’s  method  approximates  an  impulse  response 
(in  this  case  the  influence  functions)  with  a  set  of  exponential 
functions  in  a  least-squares  sense.  These  exponential  approx¬ 
imations  to  the  influence  functions  can  be  readily  transformed 
into  an  /th  order  HR  filter  of  the  form: 

/-I  //-i 

PUz)  =  2 

tmO  /  1-0 

Equating  like  powers  of  the  FIR  representation  in  Eq.  (15)  with 
the  HR  representation  of  Eq.  (16),  and  replacing  the  temporal 
argument  iT  with  subscripts,  yields  an  overdetermined  set  of 
equations  (for  R  »  /)  of  the  form: 


QUO  =  +  UOD^ 

(14b) 

(14c) 

(14d) 

Note  that  this  result  is  identical  to  that  obtained  from  full 
potential  flow  analysis  if  the  influence  functions  [HUO  and 
lUO]  are  ignored. 

Approximate  Generalized  Forces 

Many  techniques  for  approximating  the  aerodynamic  loads 
on  structures  (mainly  wing  sections)  have  been  published  in 
the  past.  Jones"  was  one  of  the  earliest,  whereas  Leishman'^ 
and  Morino  et  al.^  have  more  recently  investigated  Laplace/ 
frequency  domain  approximations  of  the  aerodynamic  loads. 
However,  the  form  of  the  aerodynamic  forces  expressed  in 
Eqs.  (12)  does  not  lend  itself  to  direct  application  of  these 
methods.  Utilization  of  the  techniques  in  Refs.  1,  9,  and  10 
would  require  transforming  the  time  domain  representations  in 
Eqs.  (12)  to  the  frequency  domain.  This  additional  computa¬ 
tional  step  would  introduce  more  error  into  the  solution.  There¬ 
fore,  a  more  direct  time  domain  approximation  technique  was 
sought. 

Since  the  generalized  aerodynamic  forces  represented  in  Eq. 
(12b)  are  functions  of  the  panel  generalized  coordinates,  they 
can  be  viewed  as  dynamic  feedback  on  the  panel.  These  feed¬ 
back  dynamics  are  characterized  by  the  influence  functions 
HUO  and  lUO,  which  represent  aerodynamic  impulse  re¬ 
sponses,  and  the  influence  coefficients  and  D^,  which  are 
instantaneous  or  feedthrough  dynamics. 

Since  there  is  no  closed-form  solution  for  the  generalized 
forces,  some  approximation  must  be  made.  The  approach  sug¬ 
gested  here  is  to  Construct  digital  filters  that  approximate  these 
dynamics.  A  finite  impulse  response  (FIR)  filter  of  the  form: 

*-1 

p{z)  =  ^hmz-‘  (15) 

f-0 

exactly  represents  a  dynamic  system  having  a  discrete  time 
impulse  response  /i(iT)  and  a  z-domain  transfer  function  P(z) 
(Ref  13).  In  this  case  the  dynamic  system  is  the  aerodynamics 


feedback  on  the  pane  iulting  in  a  complete  aeroelastic  sys¬ 
tem.  This  is  shown  schematically  (in  state-space  notation)  in 
Fig.  3  and  will  be  discussed  in  more  detail  in  the  next  section. 


Note  that  care  should  be  taken  in  the  st  ,ion  of  T.  Significant 
errors  can  result  if  the  sample  period  is  not  sufficiently  small. 
Application  of  the  bilinear  transformation  results  in  a  contin¬ 
uous  time  filter  of  the  form: 


F™(k)  =  2  diK’  / 2 


where  d,  and  are  new  coefficients. 

As  an  example,  suppose  the  values  of  the  influence  function 
are  used  to  solve  for  the  filter  coefficients  in  Eq.  (17a). 
When  provided  with  an  input  signal  that  represents  the  panel 
generalized  coordinate  this  filter  would  output  an  ap¬ 

proximate  convolution  of  q„{iT)  with  H^iiT).  Similarly,  a  fil¬ 
ter  could  be  constructed  from  and  provided  with  a  sig¬ 

nal  representing  q„{iT).  Summing  the  outputs  of  the  two  filters 
would  yield  an  approximation  to  the  dynamic  portion  of  Eq. 
(12b).  This  process  will  be  covered  in  more  detail  later. 

The  accuracy  of  this  type  of  approximation  is  demonstrated 
in  Fig.  2.  This  plot  shows  the  generalized  force  g,,  for  a  unit 
step  change  in  the  generalized  coordinate  as  represented  in 
Eq.  ( 1 2b).  The  results  are  presented  for  a  panel  of  aspect  ratio 
alb  =  1.0,  mass  ratio  /z  =  0.1,  and  M  =  2.0.  Three  curves  are 
shown,  one  representing  direct  calculation  from  the  influence 
function  while  the  others  represent  the  step  response 

of  an  8th-  and  16th-order  approximating  filter  with  a  sample 
period  of  7"  =0.1.  As  one  would  expect,  the  higher-order  filter 
yields  higher  accuracy.  Other  influence  functions  exhibit  sim¬ 
ilar  accuracy. 

An  approximate  filter  can  be  found,  as  shown  previously, 
for  each  and  A  complete  set  of  these  filters  can 

then  be  summed  according  to  Eq's.  (12a)  and  (12b).  This  would 
yield  a  dynamic  system  that  has  the  generalized  coordinates 
q„(t)  and  q„{t)  as  the  inputs,  and  the  generalized  aerodynamic 
forces  2,(0  as  the  outputs.  This  system  is  then  coupled  as 


State-Space  Formulation 

The  general  equation  of  motion  for  the  panel,  described  in 
Eq.  (3),  in  the  absence  of  structural  damping,  can  be  expressed 
in  state-space  form  as 


Xp  =  ApXp 

+  BpUp 
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92  ■  ■  ■  9n 

9i  92 

<jnV 

0 

0 

0  1  •• 

•  0~ 

0 

0 

0  0 

1 

-oA 

0 

0  0  •• 

0 

0 

3  .. 

1 

0  0  •• 

0 

_  0 

0 

-0)1  0  •• 

0_ 

-  0  0 
0  0 

0  0 

B,  =  PaVl  ^  0 

M, 


Cy  -0.02 


- —  Exact  Response 

w-  9  m  Approximate  Filter,  la  1 6 
O"  O  Approximate  Filter,  1=8 


0  0 


u,=  {Q,{t)  QiU)  •••  QMV 


1 

0  • 

••  0 

0 

0 

0 

0  • 

••  0 

1 

0 

0 

1  • 

••  0 

0 

0 

0 

0  ■■ 

•  0 

0 

1 

0 

0 

•  1 

0 

0 

0 

0  •• 

•  0 

0 

0 

Non-dimensional  time,  s 

Fig.  2  Generalized  force  step  response  for  exact  and  two  ap¬ 
proximating  filters. 
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such  that  the  output  vector,  which  is  equivalent  to  the  input  of 
the  aerodynamic  system,  is  a  reordered  version  of  the  panel 
state  vector: 

3'r={9i  9i  92  92  •••  9n  (20h) 

A  filter  of  the  form  in  Eq.  (19)  can  be  cast  in  state  variable 
form  as  follows: 


Xp  —  ApXp  -1"  SpUp 
yp  —  FpXp  -(-  GpUp 

here 

—BjBi  —By/Bi  —B,iBx 

1  0  •••  0  0 
Ap=  0  1  •  ••  0  0 


Fig.  3  State-space  schematic  of  the  complete  aeroelastic  system. 
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P" p  —  Cl\S^S  I 

Gp  =  [fli/^i] 


djB,  - 

(21e) 


The  input  is  the  appropriate  panel  generalized  coordinate  and 
the  output  is  the  approximate  convolution  of  the  generalized 
coordinate  and  the  influence  function  used  in  the  filter  calcu¬ 
lation  of  Eq.  (17a).  The  states  of  the  previous  system  are  of 
mathematical  construct  and  have  no  physical  significance. 

At  this  point  a  representation  of  Eq.  (12b)  can  be  constructed 
for  a  single  plate  mode  (m  =  n)  as  follows: 

1 )  The  appropriate  influence  functions  HppSiT)  and 

are  found  by  the  numeric  integration  of  Eqs.  (12e)  and  (12f). 

2)  Two  filters  are  constructed  by  solving  Eq.  (17a)  for  the 
filter  coefficients,  one  based  on  H,JiiT)  and  the  other  based 
on 

3)  Each  of  these  two  filters  is  transformed  to  continuous  time 
using  Eq.  (18)  and  cast  in  state  variable  form  according  to  Eq. 
(21a).  In  general,  the  state  variable  representations  associated 
with  H^{iT)  and  I^{iT)  will  be  referred  to  as  the  sets 

Bh.„  F„^,  Ch^]  and  {A,„,  F,„,  G,„],  respectively. 

4)  Finally,  the  state  variable  representations  are  combined 
according  to  Eq.  (12b)  as  follows; 

i,.  =  AaXa  +  BjUp  (22a) 

=  FpXp  +  GjUp  (22b) 


yp  =  Q™(0 


Up  = 

■ 

A„  = 

Ah^ 

*  *a 

0 

A,^J 

B.  = 

0 

F.= 

[F„_ 

F,^] 

and  the  influence  coefficients  of  Eqs.  (12c)  and  (12d)  are  in¬ 
corporated  as  follows: 


Go  - 


+  Omm] 


The  previous  representation  would  be  sufficient  if  the  panel 
model  had  only  one  mode.  This  would  result  in  only  one  term 
in  the  sum  given  by  Eq.  (12a).  Multiple  panel  modes  require 
the  extension  of  the  previous  representation  to  include  all  ap¬ 
propriate  influence  functions.  This  more  general  representation 
is  as  follows; 

Xa  =  AaXp  +  (23a) 

yo  =  FpXp  +  Gju.^  (23b) 
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and  the  influence  coefficients  are  incorporated  as  follows: 

G'h„  =  G„^  +  S^  .(23i) 

GL  =  +  Dpp,  (23j) 

The  panel  and  aerodynamic  state-space  models  can  then  be 
coupled  together  as  shown  in  Fig.  3. 

TTie  size  of  the  state-space  model  can  be  further  reduced  by 
taking  advantage  of  certain  symmetries  in  the  influence  func¬ 
tions.  The  greatest  reduction  can  be  had  by  noting  that  = 
and  Furthermore,  if  more  than  one  expansion 

mode  is  used  in  the  spanwise  direction  for  a  simply  supported 
plate,  the  influence  functions  can  be  indexed  such  that  W™,  — > 
Hp^r,  and  7™  — >  In  this  case,  those  influence  functions 

where  —  sj  is  odd,  are  zero. 

Results 

To  validate  the  coupled  model,  the  system  parameters  were 
varied  and  the  flutter  boundaries  were  compared  to  results 
from  the  literature.  Flutter  boundaries  were  found  by  increas¬ 
ing  the  nondimensional  dynamic  pressure  A  and  calculating  the 
coupled  system  eigenvalues.  Flutter  instability  occurs  when  at 
least  one  of  the  system  eigenvalues  lies  in  the  right  half-plane. 

This  is  demonstrated  in  Fig.  4,  which  shows  the  distribution 
of  the  system  eigenvalues.  Note  the  distribution  of  the  ap¬ 
proximate  aerodynamic  eigenvalues.  Those  eigenvalues  asso¬ 
ciated  with  the  panel  modes  are  contained  in  the  zoom  region 
indicated  in  Fig.  4.  These  results  are  for  a  panel  with  db  = 
1.0,  /j.  =  0.1,  Af  =  2.0,  and  using  a  four-mode  expansion  [i.e., 
N  =  4  in  Eq.  (2)].  In  this  case,  the  frequency  is  nondimen- 
sionalized  such  that  k  -  (tocz/(/„)(A//x)''^. 

Figure  5  is  an  expanded  view  of  the  zoom  region  of  Fig.  4 
and  depicts  the  pole  migration  with  increasing  dynamic  pres¬ 
sure.  In  Fig.  5  eigenvalues  at  different  dynamic  pressure  are 
connected  with  a  line.  Note  the  path  of  the  complex  conjugate 
eigenvalue  pair  as  they  enter  the  right  half-plane.  Also  note 
that  the  only  eigenvalues  that  experience  significant  migration 
are  the  complex  conjugate  pairs  that  are  associated  with  the 
panel  modes.  One  of  these  two  migrating  pairs  is  the  eigen¬ 
value  that  becomes  unstable. 

Another  important  piece  of  information  that  can  be  obtained 
from  the  eigenvalue  locus  in  Fig.  5  is  the  frequency  of  flutter. 
In  Fig.  5  the  frequency  of  flutter  is  between  the  first  and  second 
natural  frequencies  of  the  panel.  This  is  consistent  with  the 
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Fig.  4  Locus  of  system  eigenvalues  with  varying  dynamic 
pressure. 
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Fig.  S  Close-up  locus  of  system  eigenvalues  with  varying  dy¬ 
namic  pressure,  Af  =  2.0. 


Mach  Number 

Fig.  6  Flutter  dynamic  pressure  vs  Mach  number  for  a/b  =  0. 

A  comparison  of  flutter  dynamic  pressure  vs  Mach  number 
is  shown  in  Fig.  6.  Results  are  illustrated  for  a  panel  with 
/i  =  0. 1  and  alb  =  0  (i.e.,  an  infinitely  wide  panel).  Also  in¬ 
cluded  in  the  figure  are  the  flutter  boundaries  for  a  piston  the¬ 
ory  model  and  for  an  approximate,  quasisteady  aerodynamic 
flutter  boundary  described  by  Hedgepeth.’ 
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Fig.  7  Flutter  dynamic  pressure  vs  Mach  number  for  atb  =  1. 
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Fig.  8  Flutter  dynamic  pressure  vs  damping  ratio. 

Note  in  Fig.  6  that  the  three  theories  approach  each  other 
asymptotically  as  the  Mach  number  increases.  However,  for 
Mach  numbers  below  1.5,  the  theory  incorporating  linearized 
potential  flow  aerodynamics  predicts  flutter  boundaries  well 
below  the  others.  Full  potential  flow  theory  predicts  a  mini¬ 
mum  flutter  dynamic  pressure  of  zero  at  a  Mach  number  of 
1.0.  These  results  are  consistent  with  those  obtained  by 
Dowell.*’ 

A  similar  flutter  boundary,  for  a  square  panel  (i.e.,  atb  =  1 .0) 
with  ^  =  0.1  and  using  a  four-mode  expansion,  is  shown  in 
Fig.  7.  These  results  are  also  consistent  with  those  obtained 
by  Dowell.  Note  in  this  case  the  minimum  flutter  dynamic 
pressure  is  nonzero. 

A  final  and  interesting  result  is  the  effect  of  proportional 
modal  damping  on  the  flutter  dynamic  pressure.  Figure  8 
shows  the  flutter  dynamic  pressure,  normalized  with  respect  to 
the  undamped  flutter  dynamic  pressure,  as  a  function  of  the 
modal  damping  ratio  for  three  different  Mach  numbers.  Pro¬ 
portional  damping  is  modeled  as  an  additional  term  of  the  form 
to  Eq.  (3).  Note  that  proportional  damping  has  the 
greatest  effect  on  flutter  dynamic  pressure  at  lower  Mach 
numbers. 

Conclusions 

A  new  method  for  modeling  the  aerodynamic  loading  of 
panels  that  incorporates  linearized  potential  flow  aerodynamics 
has  been  developed.  This  was  accomplished  by  approximating 
the  aerodynamic  generalized  forces  on  the  panel  with  discrete 
time  filters.  These  filters  were  coupled  as  feedback  to  the  in 
vacuo  panel  model.  The  accuracy  of  the  method  was  demon- 
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strated  by  comparing  results  of  this  model  with  results  ob¬ 
tained  in  previous  work.  Future  work  will  use  this  new  aero- 
elastic  model  in  optimal  control  studies. 
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